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Abstract. Imaging the magnetic fields around a non-magnetic impurity can 
provide a clear benchmark for quantifying the degree of magnetic frustration. 
Focusing on the strongly frustrated J1-J2 model and the spatially anisotropic 
Jia-J\b~J2 model, very distinct low energy behaviors reflect different levels of 
magnetic frustration. In the J1-J2 model, bound magnons appear trapped near 
the impurity in the ground state and strongly reduce the ordered moments for 
sites proximal to the impurity. In contrast, local moments in the J\a-J\b-Jl model 
are enhanced on the impurity neighboring sites. These theoretical predictions 
can be probed by experiments such as nuclear magnetic resonance and scanning 
tunneling microscopy, and the results can elucidate the role of frustration in 
antiferromagncts and help narrow the possible models to understand magnetism 
in the iron pnictdies. 
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One of the major thrusts in condensed matter physics concerns the interplay of 
strong frustration and quantum fluctuations in magnetic systems pQ. These systems 
can exhibit (a) novel phases of matter lacking conventional order, (b) a high degree of 
residual entropy at low temperatures with local constraints leading to emergent gauge 
fields [2] , and (c) new types of particles with unusual quantum numbers and dynamics 
not captured within the established framework of quasi-particles [3j. But, how can 
such magnetic frustration be recognized and quantified experimentally? 

The ratio of the Curie- Weiss temperature to the magnetic ordering temperature 
has served as an important metric [1], though it cannot distinguish between the role 
of fluctuations due to reduced dimensionality and that due to frustration. The spin- 
wave spectrum can at times lead uniquely to the underlying magnetic interactions, 
though there are cases where this also can be ambiguous. In some classical systems, 
such as spin- ice, the residual entropy can be measured directly by specific heat [5]. 
Various ratios associated with the peak in the uniform susceptibility also can quantify 
frustration [6]. However, thermodynamic measurements are not easy to interpret when 
other low energy degrees of freedom, besides spins, are present. Since frustration is a 
local property, measurements which can infer it locally would be useful. 

Studying antiferromagnetism, magnetic frustration, and quantum criticality 
in the iron pnictides is crucial for a comprehensive understanding of their 
superconductivity. This also can lead to the development of a broader framework 
for unconventional superconductivity in materials such as cuprates, ruthenates, heavy 
fermions, and organic superconductors [7]. In this letter, we explore how imaging the 
magnetic fields around a non-magnetic impurity can provide a clear benchmark for 
quantifying the degree of magnetic frustration. While the pnictides provide the context 
for this study, it should have a broader applicability in frustrated spin systems [8]. 

Theories for the collinear antiferromagnetic (AF) phase in iron pnictides vary 
greatly in the importance of weak versus strong coupling [9] [TOj [TTJ [12j [13], and 
in the amount of magnetic frustration and proximity to a quantum critical point 
(QCP) p~4] [15]. There are various reasons why a local moment perspective is still 
relevant in these systems. Parametrization of the magnetic behavior by effective local 
moment models, even for metals, should remain valid as long as there are long-lived 
magnon-like excitations |161 117) . Most importantly, however, is that local moment 
models are still being widely used in experiments to explain magnetic properties in 
these materials. These include the spatially anisotropic J\ a -J\b~J2 model [l"8l IPO] [20] . 
where the coupling is ferromagnetic (FM) in one direction and AF in the other, and 
the strongly frustrated J\-Ji model [2T] [33] [23] , where collinear AF order arises as a 
result of a subtle selection via order by disorder [24] . 

Neutron scattering on CaFe2As2 indicates that the spin wave energy is a 
maximum at momentum transfer (tt, tt) [20] , which favors the J\ a -J\b~J2 scenario. 
However, the results have been disputed [25] and alternative experiments would be 
useful. A particularly good example of competing models is the iron-chalcogcnidc 
superconductors, Fei+yTei-^Sea;, where the (tt/2,tt/2) AF order can be obtained 
by either invoking strong frustration in a J1-J2-J3 model [26] or utilizing a model 
with strong spatial anisotropy arising from orbital order [27] . Both models can 
lead to virtually indistinguishable spin-wave spectra. Hence, further experiments are 
necessary to clarify the nature of the magnetic phase. 

Leaving aside the question of where such strong in-plane anisotropy can come 
from [28l [29l l30l l3il I3"2l [33] . one can ask whether experiments can directly sense 
the degree of frustration in the system. In the pnictides, one way to distinguish 
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Figure 1. (Color online) (a) Schematic for the collinear AF state on a 2D square 
lattice, (b) The position labelings for selected sites. The impurity is placed at 
(0,0). 



models with different levels of frustration is to study the response to non-magnetic 
impurities [53], such as those caused by Zn substitution [35 . Such impurities can 
suppress or enhance antiferromagnetism in their neighborhood, and lead to bound 
magnon states with different degrees of anisotropic perturbations. 

Using a T- matrix approach within linear spin- wave theory |36j . in the present 
study we find that the static and dynamic perturbation of collinear antiferromagnetism 
near a non-magnetic impurity are markedly different in the two models. In the 
frustrated J\- J% model, a non-magnetic impurity strongly reduces its neighboring local 
moments and overturned dynamical spins appear close to zero energy, rendering non- 
trivial physics in proximity to the QCP. In contrast, the Jia-Ju-Ji model produces 
enhanced local moments on the impurity neighboring sites, a behavior similar to the 
Heisenberg model [3B]. In both models, the disturbance on the magnon local density 
of states (LDOS) heals quickly along the AF direction but is extended in the FM 
direction. This leads to spin fluctuations with an anisotropic stripe pattern. 

We consider a spin Hamiltonian on a 2D square lattice 

Hq = ( Ji a S r ■ S r +£ + JibS. r 



-Si ■ Sj 



(1) 



where ...» represents a second nearest-neighbor sum. When Jit, = J\ a the 
Hamiltonian reduces to the J\-Ji model. For a collinear AF order, we introduce 
the real space Holstein-Primakoff bosons in the linear spin wave approximation: 



5' 
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2Sa l: and 5? 
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2Sb\. We use the notation 



3 ' 3 ' 

that i belongs to sub-lattice A (spin up), and j belongs to sub-lattice B (spin down) 
[Fig. 1(a)]. A standard Bogoliubov-de Gennes transformation diagonalizes Hq with 
eigen-energies E k = J 2 SZy/A 2 k - Bl, where k is defined over the reduced (magnetic) 
Brillouin zone. Here Ak = 1 + a + (3(cosk y — 1), and B^ = cosfc x (cos/c a -I- a). Z is 
the coordination number, and a = Ji a /2J2, f) = Ju/2J2- 

For the J1-J2 model, we focus on J2 = J\ 37 . In the collinear AF state, the 
results for various J2/J1 ratios are qualitatively the same, except for J2/J1 very close 
to the QCP as we discuss below. For the J\ a -J\b-J2 model, we use Ju, = — 0.1Ji a 
and J 2 = 0.4J la , appropriate for CaFe 2 As 2 [10] ■ In the latter case the system is 
unfrustrated even if (Ji Q + Ju) is comparable to 2J 2 . 

The effects of a single non-magnetic impurity are introduced by removing the 
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Figure 2. (Color online) J1-J2 model: The imaginary parts of the Green's 
functions related to the magnon LDOS at selected positions. The magnetic 
disturbance heals quickly along the AF direction but extends more in the FM 
direction. 



exchange interactions between the impurity and its first and second nearest neighbor- 
sites. The impurity problem is then solved conveniently by defining the Green's 
function matrix: 

r ~ (f) = ( -i(Tal(t) ail (0)) -i{Tat(t)bj, (Q)) \ 
r ' r ' {) \ -ipbjWavW) -iiTb^b),^)) ) ■ [) 

The bare Green's functions in momentum and frequency space are 



TQJ + J 2 SZA k 

(22) v ' n ' n JJ-1.0+ to 2 — Ef. + in 



(G°k, k ') 11 M = h,k' lim 



- (") = h,k> Hm 2 Jz ^f f . ■ (3) 
(21) 77^0+ oj z — hf, + 17] 

The dressed Green's functions G in the presence of an impurity can be calculated 
exactly in the thermodynamic limit (within spin wave approximation) by the T-matrix 
formalism 36 . This theoretical treatment has ensured that free bosons associated with 
the impurity site never hop and thus avoids spurious zero modes |38j . The physical 
quantities of interest include the magnetic moments, as well as the magnon LDOS, 
which can be obtained from the imaginary part of the Green's function [36] . 

One of our most interesting findings is the presence of overturned dynamical spins 
near the impurity in the ground state for the J\-Ji model, but not for the J\ a -J\b~J2 
models. In the J%- J2 model, the magnon LDOS for spin sub-lattice B [left panels of 
Fig. 2] shows more low energy features, and localized magnon excitations appear near 
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Figure 3. (Color online) J\a-Jib~J2 model: The imaginary parts of the Green's 
functions related to the magnon LDOS at selected positions. The magnetic 
disturbance heals quickly along the AF direction but extends more in the FM 
direction. 



the impurity at finite, positive energies 0. In contrast, the LDOS in the J\ a -J\b-J2 
model shows features only at relatively high energies close to the van Hove singularity 
of the homogeneous host system, as shown in Fig. [3] In both models, the disturbance 
on the magnon LDOS heals quickly along the AF direction but extends more in the 
FM direction. 

We also find a strong reduction in the magnetic moments near the impurity sites 
in the frustrated model, but not in the unfrustrated case. Figure 0] shows the change 
in magnitude of the magnetic moments (|(>S£)| — \(S^} \)/\(S^} \ for the two models. 
In the frustrated J1-J2 model, at site r — (1,0) (belonging to a different spin sub- 
lattice from the impurity), the magnitude of the moment strongly decreases, indicating 
enhanced quantum fluctuations [Fig. 4(a)]. On the other hand, in the unfrustrated 
J\a-J\b-J2 model the magnitude of the moment increases [Fig. 4(b)], a behavior more 
similar to the Heisenberg model [36]. Moreover, while the change in magnitude of the 
moments in the J\ a -Jib-J2 model is ~ 10%, this change in the J1-J2 model is more 
substantial. In both models, the change of the moments is damped quickly away from 
the impurity by a few lattice sites, but exhibits spin fluctuations of an anisotropic 
stripe pattern. Density functional calculations on cobalt doped pnictides also found 
a similar behavior [5T]. We emphasize that these fluctuations are associated with the 

I There are artificial zero modes at k = (it, it) in the linear spin wave spectra in the J1-J2 model, 
which is known to be gapped by quantum fluctuations or spin anisotropy 39 , 40 . By making J\ a and 
Jib unequal, we have verified that our results change only slightly when a spin-wave gap at k = (n, w) 
is gradually turned on. 
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Figure 4. (Color online) Changes in the magnetic moment for (a) J1-J2 
model and (b) J\ a -Jib~J2 model. A spline interpolation is used for both plots. 
The moment at site r = (1,0) [{S^ , 1Q ,)| is strongly reduced in the J1-J2 
model indicating a stronger quantum fluctuation, while for the Jia-Jif,- J2 model 
\(S*_q j)| is slightly enhanced. The change in magnitude of the moments in the 
Jla-J\b~J2 model is ~ 10%, but this change in the J1-J2 model is more substantial. 
Both models exhibit spin fluctuations with an anisotropic stripe pattern. 



spin variables, which must be distinguished from the low-energy charge excitations in 
metals. 

Aside from changing the magnon LDOS and the magnetic moments, non-magnetic 
impurities can also result in an instability of the ground state 42 . In the spin 
wave calculation, to sustain a collinear AF ground state in the J\-J% model requires 
J2 > 0.5 Ji. With the presence of a non-magnetic impurity, our calculations for 
J2/J1 < 0.6 indicate that the moments are drastically reduced, and the changes in 
moments exceed the bare values. In this case, the system no longer favors a collinear 
AF state and the spin wave calculation breaks down. In the pnictides, the J1-J2 
model requires fine-tuning to the QCP to explain the reduced ordered moment [39J. 
Therefore, a competing phase other than the collinear AF state may result upon 
doping by non-magnetic impurities. 

In principle, the ordered local moments can be measured directly in nuclear 
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magnetic resonance (NMR), where they result in a local magnetic field on the nucleus 
and cause a shift in its resonance frequency [53] • Our calculations for non-magnetic 
impurities (such as Zn substituting Fe) indicate that the J\-J% model produces a more 
substantial broadening of the magnetic field distribution, while this inhomogeneous 
broadening in the Jia-Jib-Ji model is much smaller. On the other hand, in both cases 
the presence of line defects such as those coming from domain boundaries would tend 
to spread the distribution. These issues as well as confirmation of the anisotropic 
disturbance of the collinear AF state deserve further experimental attention. We last 
note that recent NMR measurements on Ni doped BaFe2As2 found that the system 
remains in a collinear AF state upon doping [44] . The experiments also show a 
considerable broadening of the field distribution. However, it is not clear to what 
extent Ni or Co dopants can be treated as localized non- magnetic impurities [351 HB] . 

Spin-resolved scanning tunneling microscopy (STM) can potentially measure 
the full frequency dependence of the spin structure [?71 HE]- One complication is 
that there are other low energy (charge) degrees of freedom present besides spin. 
Indeed, STM measurements for low energy quasi-particles on the pnictides show 
unidirectional nanostructures weakly pinned by dopant atoms, strongly suggestive 
of broken tetragonal symmetry and in-plane anisotropy [49[ [50], On the other hand, 
the spin-resolved magnetic properties may be embedded under signals from the charge 
degrees of freedom, and proper care is needed to separate the tunneling conductance 
between localized spins and itinerant charges. A spin-polarized STM with atomic 
resolution can filter out the magnon-like excitations, whose projection is fully spin 
polarized at every site, from those itinerant charge degrees of freedom, which may be 
unpolarized or weakly spin polarized. 

In summary, we have studied the responses to a non-magnetic impurity for the 
J1-J2 and J\a-J\b-J2 models in the collinear AF state. We found that frustration 
in the J1-J2 model results in bound magnons at very low energies with concomitant 
reduction of the magnetic moments of the impurity neighboring sites. In contrast, the 
J\a-J\b-J2 model shows behavior more similar to the square-lattice Heisenberg model, 
where local moments are enhanced on the sites proximal to the impurity. We made 
clear experimental predictions for magnon LDOS and real space spin textures that 
potentially can be measured by NMR, STM or other resonance probes. The results 
can benchmark strong frustration and proximity to a QCP in the system under study. 

We last note that the results discussed above should qualitatively carry over 
to other spin systems with different AF order, such as the (ir/2,ir/2) order in the 
Fei+yTei-^Se^ materials. Due to differences of frustration intrinsic to the proposed 
models 27, 28], the magnon LDOS and the response of the AF hosts to non-magnetic 
impurities can be rather different. Calculations related to this interesting extension 
are an area for future work. 
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